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> ; 1 INTRODUCTION 

I In the classical theory of structure formation, it was gener- 
$_( ' ally believed that the gravity is fully responsible for the for- 
, mation and evolution of the large scale structure in the uni- 
verse. Recently, however, it has been realized that the overall 
characteristics of the large scale structures cannot be un- 
derstood only in terms of the gravitational influence. It has 
been pointed out that the tidal field provides a driving force 
in establishing the o bserved large sca l e filamentary struc- 
tures in the universe jBond et aU h99&. 'L ee fc Evrardll2007l: 
Park fc Led l2007bl : iHahn et al.r 2007; A ragon-Calvo et al l 
cosmic voids are most vulnerable to the tidal 



ABSTRACT 

The size distribution of mini-filaments in voids has been derived from the Millennium 
Run halo catalogs at redshifts z = 0, 0.5, 1 and 2. It is assumed that the primordial 
tidal field originated the presence of filamentary substructures in voids and that the 
void filaments have evolved only little, keeping the initial memory of the primordial 
tidal field. Applying the filament-finding algorithm based on the minimal spanning 
tree (MST) technique to the Millennium voids, we identify the mini-filaments running 
through voids and measure their sizes at each redshift. Then, we calculate the comoving 
number density of void filaments as a function of their sizes in the logarithmic interval 
and determine an analytic fitting function for it. It is found that the size distribution 
of void mini- filaments in the logarithmic interval, dN/d\ogS, has an almost universal 
shape, insensitive to the redshift: In the short-size section it is well approximated as 
a power-law, dN/dlogS ~ S, while in the long-size section it decreases exponentially 
as dN/dlogS ~ exp(— S*"). We expect that the universal size distribution of void 
filaments may provide a useful cosniological probe without resorting to the rms density 
fluctuations. 

Key words: cosmology:theory — large-scale structure of universe 
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influence from the surrounding matter distribution due 
to their extreme low-dens i tv llSahni fc ShandarinI [l99^: 
Shandarin et~aLll2004 l2006l : iLee fc Parkll2006l : iPark fc Led 



2007bl ). 



The tidal squeezing and distortion effect tends to 
deviate the void shapes from spherical symmetry and could 
sometimes le ad even to the c o llapse and disappearance 
of the voids llSahni et all Il994l : ISahni fc ShandarinI 1 19961 : 
IShandarin et al.ll2004l ')! 

The presence of the filamentary structures in voids 
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marks the most striking evidence for the strong tidal in- 
fluence on voids. The anisotropy in the spatial distribution 
of void halos is induced by th eir alignments with the princi- 
pal axes of the tidal tensors (|Peeblesll200ll ). Since the void 
halos ev olve very little and most particles remain primordial 
in voids (lEinastoll2006l ), the void mini-filaments are pristine, 
keeping well the initial memory of the spatial coherence of 
the primordial tidal effect, unlike the large scale filaments 
which also originated from the large-scale coherence of the 
primordial tidal field but have undergone highly nonlinear 
processes in the subsequent evolution. Hence, the void fila- 
ments are the most idealistic probe of the primordial tidal 
field and its effect on voids. 

The size distribution of void filaments should be a func- 
tion of the strength of the tidal effect and the host void sizes. 
The maximum size of a void filament cannot exceed the size 
of its host void. But, even when a host void has a large size, 
it would not have long-size filaments if there is no tidal in- 
fluence. Here, we attempt to derive the size distributions of 
the void filaments at various redshifts and to explore their 
statistical properties. Before presenting our results, however, 
we would like to caution the readers for the difficulty in us- 
ing the void filaments as a probe of the primordial tidal 
field. Unlike the galaxy clusters, there is no unique way to 
define voids and their filaments. Recently, the seminal work 
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Figure 1. The statistical significance P{sc) vs. the minimum 
void-size threshold Sc at 2; = 0, 0.5, 1 and 2. 

of IColberg et alj l|2008h compared thoroughly various void- 
finding algorithms and showed clearly that different void 
finders different void characteristics, under-density profiles, 
void galaxies etc. Thus, our results are likely to be dependent 
on our specific choice of the void-finder. The organization of 
this paper is as follows. In ^ we briefiy describe the void 
catalog was obtained from the Millennium Run simulations. 
In we explain how the void filaments are identified from 
the void catalog. In i|4l we derive numerically the size distri- 
butions of void filaments and determine an analytic fitting 
formula for it. In SJSI we discuss the final results and assess 
the future work. 



2 OVERVIEW OF THE VOID 
IDENTIFICATION 

In our previous work l|Park fc "Le3l2007bD . we have obtained 
a sample of voids at four different redshifts, z = 0, 0.5, 1 
and 2, from the Millennium Run Simulation (|Springel et al.l 
I2OO5I ) for a flat ACDM cosmology with the key cosmological 
parameters Q.„, = 0.25, Hi, = 0.045, h = 0.73, = 0.75, 
n = 1, and cts = 0.9, which we use throughout this paper. 
Here, we briefiy summarize the void-identification process 
as an overview. 

We first selected those halos in the Millennium data 
which consist of more than 50 particles since those halos 
composed of less than 50 halos are statistically unreliable 
(in private communication with V. Springel). Beside, it is 
interestingly found that when the particle number thresh- 
old of Nc = 50 is used, the statistical properties of the halo 
voids are quite similar to that of the galaxy voids from the 
Millennium Sloan Digital Sky Survey (SDSS) mock cata- 
log with the magnitude completene ss limit. We applied the 
void-finding algorithm developed bv lHovle fc VogelevI l|2002l . 
hereafter HV02) to the sample of the selected halos at each 
redshift, separately. 

To apply the HV02 algorithm, it is required first to 
set the values of two parameters: the wall/field criterion 
I and the minimum void-size threshold Sc. The value of 
the wall/field criterion was calculated as I = + 3ct3/2, 



Table 1 . The total number of the Millennium voids , the mean 
density contrast Sy, and the mean effective radius Ru a,t z = 
0,0.5, 1 and 2. 



z 




Sv 


Rv 





29186 


-0.90 


9.51 


0.5 


27652 


-0.90 


9.42 


1 


26615 


-0.89 


9.39 


2 


27951 


-0.89 


9.42 



where da is the mean distance to the third nearest neigh- 
bor halo and is its standard deviation. If is found that 
/ = 2.78, 2.70, 2.67, and 2.75/i"^Mpc, at z = 0,0.5,1 and 
2, respectively. To determi ne the size threshold , we tested 
the statistical significance ||E1-Ad fc Piranl l 19971) that is de- 
fined as P{sc) = 1 — A''gap(sc)/A''void (sc) where A'"„ow(sc) 
and Ngap{sc) are the numbers of the voids found in the real 
sample and in 10 independent random samples, respectively. 
This test basically allows us to discriminate the true voids 
from the Poisson gaps. Fig. [1] plots -P(sc) as a function of 
the size threshold, Sc at z = 0, 0.5, 1 and 2 (solid, dashed, 
dot-dashed and dashed line, respectively) . Noting that when 
Sc = 5.5/i~^Mpc, -P(sc) reaches 0.98 at all redshifts, we set 
the void size threshold at 5.5/i~^Mpc. 

Table [T] lists the statistical properties of the voids at 
each redshift. It is worth mentioning here that these proper- 
ties will depend on t he void-finding algorithm. According to 
IColberg et al.l (|2008l ). the HV02 algorithm appears to iden- 
tify as the void a much larger region (and thus identify many 
more galaxies as void-galaxies) than most of the other algo- 
rithms. Thus, the mean effective comoving radius Rv and the 
number of void galaxies A'^„ listed in Table [1] are expected to 
be larger than for the cases of most of the other algorithms. 
It also indicates that the number of void filaments and their 
sizes would be larger when the HV02 algorithm is used. 



3 THE VOID FILAMENTS IN A ACDM 
UNIVERSE 

We em ploy the filament-finding algorithm of iBarrow et al.l 
(| 19851 ) that is based on the minimal spanning tree (MST) 
scheme to extract the intrinsic linear patterns from the spa- 
tial distribution of the halos in the selected large voids which 
contain more than 30 halos. There are three key concepts for 
the MST scheme: node, edge and fc-branch. A node means 
the point (i.e., halo), an edge means a straight line that con- 
nects two nodes, and a fc-branch means a path consisting of 
fc-edges with a leaf node only at one end. That is, one end 
of a fc-branch is free while the other end is an intersection of 
the edges. The MST of A'^ data points represents a unique 
network connected by A'^ — 1 edges of minimum total length 
without containing a closed loop. 

To construct MSTs from the spatial distribution of the 
void halos, we follow the simplest Prim process: First, an ar- 
bitrary node (halo) is chosen as a starting point. Then the 
starting point is connected to its nearest neighbor along a 
straight line (edge). The edge composed of two nodes repre- 
sents the first partial tree. The next nearest node is deter- 
mined and added to the partial tree with an edge, via which 
the partial tree is extended. After all nodes are connected to 
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the tree according to the above prescription, the MST con- 
struction is completed. The dominant linear patterns, i.e. 
filaments, are now extracted fr om the MSTs by ta king the 
following MST-reduction steps (jBarrow et al.|[l985l ): 

(i) Pruning: An MST is pruned to the p level by elimi- 
nating all the fc-branches with k^p. 

(ii) Separating: All the edges longer than a given cut-off 
length Ic are removed. 

The Pruning process allows us to find the main stems, re- 
moving the minor twigs that hardly contribute to the main 
structural patterns. In other words, the superfluous small- 
scale noises from the MSTs are minimized and the promi- 
nent filamentary patterns are highlighted by the Pruning 
process. On the other hand, the Separating process breaks 
apart the MSTs into distinct pieces, cutting off unphysical 
long linkages. Both the pruning level p and the cut-off length 
Ic quantify the degree of the MST-reduction. 



In most of the previous works where the large- 
scale filaments were the targets to find by means of 
the MST technique, the values of p and Ic were empir- 



icallv determined as d = 9 or 10 (IBarrow et al. 


19851: 


Bhavsar & Lina 


1988al bf) and L = fl (iBarrow et all 


1983: 


Bhavsar & Lina 


1988bl:lPlionis et al.lll992l:|Pearson & Coles 


19951: Krzewina & Saslawl 19961: Coles et al.l 1998h where I 



is the mean edge-length of the unreduced tree. And there 
were some authors who preferred Ic = I + 2ai or 1^ = I + cti 
l|Bhavsar fc Lingj|l988al : lzucca et al.lll99ll '). Unlike the pre- 
vious works, however, our target is not the large-scale fila- 
ments but the mini-filaments in void regions. Accordingly, 
the relevant values of p and Ic are supposed to be different in 
our case. As a practical strategy to determine the values of 
p and Ic for the void filaments, we investigate how the final 
result (i.e., the size distribution of void filaments) changes 
with p and Ic and look for those values at which the final 
result stabilizes. 

With the help of this strategy, we set p — 4 and 
Ic = I -\- ai for the void filaments. Here, the range of p is 
obtained empirically from the distribution of k values of 
the unreduced trees. By setting p = 4, we use only those 



Table 2. The total number of voids with more than 30 hales (Nv), 
the total number of void filaments (Nf), the range of the filament 
numbers per a void (nmax), the mean number of filaments per 
a void (fif), the mean length of void filaments (Rf) in unit of 
h~^Mpc, and the mean linearity of void filaments (Rl) at 2 = 
0, 0.5, 1, and 2. 



2 






nf 


nf 


[/i-iMpc] 


Rl 





5604 


22697 


[1,22] 


4 


15.06 


0.72 


0.5 


5992 


24130 


[1,16] 


4 


14.79 


0.72 


1 


6039 


24095 


[1,16] 


4 


14.69 


0.72 


2 


5349 


20514 


[1,16] 


4 


15.11 


0.71 



branches which consist of at least 5 edges as the main stems 
among the fc-branches. The detailed justification of p = 4 
is given in ^ As for Ic = I + ai, we examine plenty of 
the Millennium voids and find that the Millennium void ha- 
los have maximum number of branches when Ic — I + ui 
ICraham et al.lll995l : iBastian et al.ll2007h . After completing 
the MST-reduction process, we end up with a final sample of 
the void filaments (reduced branches) with more than one 
edges. Fig. [2] depicts the spatial distribution of the mini- 
filaments in a Millennium void on a two dimensional pro- 
jected plane. Table [5] lists the statistical properties of the 
void-filaments. Here, the linearity of a filament, Rl, is de- 
fined by the ratio of the straight distance between the two 
end nodes of the filament (end-to-end separation) to the 
sum of the len gths of all the edges constituting the filament 
( Ba rrow et al.l 119851. Basically, it represents the degree of 
the straightness of a filament: The closer to unity the value 
of Rl is, the more straight a given filament is. 

Figure [3] plots the number distribution of void fila- 
ments, Nf{RL;z) vs. the linearity Rl at z = 0, 0.5, 1 
and 2 as solid, dashed, dot-dashed, and dotted histogram, 
respectively, with single-edge filaments included (left) and 
excluded (right). As can be seen, when the single-edge fila- 
ments are included, the distribution has a peak around unity. 
In contrast, when the single-edge filaments are excluded, the 
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Figure 3. Number distribution of void filaments as a function of the linearity, i?^ at 2; = 0, 0.5, 1, and 2 as solid, dashed, dot-dashed, 
and dotted line, respectively, with single edges included (loft) and excluded (right). 




Figure 4. Number distribution of voids as a function of the num- 
ber of their member filaments at z = 0, 0.5, 1, and 2 as solid, 
dashed, dot-dashed, and dotted line, respectively. 



peak occurs at Rl = 0.63 and the mean linearity is reduced 
to Rl ~ 0.64. This indicates that except for the single-edge 
filaments most of the void filaments are not straight but 
curved. It can be readily understood since the void filaments 
are usually located in the boundary of voids. 

Figure [3] plots the number distribution of voids, 
Nv{Nf; z), as a function of the number of member filaments, 
Nf, at 2 = 0, 0.5, 1 and 2 as solid, dashed, dot-dashed, and 
dotted histogram, respectively. Most of the voids have less 
than 10 member filaments. 



THE SIZE DISTRIBUTION OF VOID 
MINI-FILAMENTS 



We first adopt the definition of lColberd l|2007l ) for the fila- 
ment size, S, and measure S of a given filament as the length 
of the diagonal line of a cuboid that exactly fits the spatial 



extent of the filament nodes. Let (xmin, Xmax), (ymin, i/max), 
and (zmin,Zmax) represent the full ranges of the x, y and 
z positions of the filament nodes, respectively. Then, S is 
determined as 



S = 



{yn 



•(1) 



By equation ((TJ, we measure the size of every void filament 
and determine the number density of void filaments as a 
function of log 5*. Fig. [5] plots the size distribution of void 
filaments, dN/dlog S, for the six different cases of the prun- 
ing level p at z — 0, 0.5, 1, and 2 (solid, dashed, dot-dashed, 
and dotted line, respectively). When p increases from 1 to 
3, there are noticeable changes in the shape of dN/dlogS. 
Whereas, when p increases from 4 to 6, there is only little 
change in the shape of dN/dlogS. It indicates that the on- 
set of the stabilization of dN/d\ogS occurs at p = 4. No 
branches are further pruned and thus the main structures 
of the MSTs are retained after p = 4. This trend justifies 
our choice of the pruning level p = 4 in §3. As can be seen 
in Fig. O the distribution dN/d\ogS appears insensitive to 
z. The four lines in each panel behave very similarly, which 
implies that it may be possible to express dN/d\ogS in a 
universal form. Given the fact that the sizes of the void fil- 
aments depend on the sizes of their host voids, we use the 
rescaled filament size, v = S/Rv where Rv is the effective 
comoving radius of a host void that is calculated by means 
of the M onte Carlo integrati on method as described in HV02 
(see also iPark fcTe3l2007bl ) . 

It is found that the rescaled size distribution of void 
filaments dN/ d log f at all redshifts can be well fitted to the 
following analytic formula: 



f{u) = Ai^exp i^ — Bu^ 



(2) 



where A and B are the fitting parameters whose best-fit val- 
ues are found through the ^'^-minimization. Fig. |6] plots the 
rescaled size distributions of void filaments at four redshifts 
and compares them with the analytic fitting formula with 
the best-fit parameters. As can be seen, the numerical re- 
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Figure 5. Number density of void-filaments per unit volume in a logarithmic size bin for the six different cases of the pruning level, 
at z = 0, 0.5, 1 and 2 (solid, dashed, dot-dashed, and dotted line, respectively). 
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Figure 6. Number density of void filaments as a function of the 
rescaled size at z = 0, 0.5, 1 and 2 (open circles, squares, tri- 
angles, and diamonds, respectively). The solid line represents a 
fitting function (eq. [2]) with the best-fit values of A and B. The 
filament size is measured as the spatial extent of the filament 
nodes. 



suits at four redshifts are indeed in good agreement with the 
fitting model. 

Yet, a careful reader may well suspect that the univer- 
sality and the functional form of dN/ d log v may depend on 
the definition of the filament size. To test how the fitting 
model and the behavior of dN/dlog v changes with the def- 
inition of the filament size, we redefine the size of a filament 
as its total length, St, (the sum of the lengths of all edges 
that constitute the filament) and repeat the whole process. 
It is found that dN/ d log f is still almost universal but the 
functional form of the fitting model is different. 

f{u) = Auex.p{-Bv) . (3) 

Fig. [7] plots the same as in Fig. [6] but for the case that the 
filament size is measured as the total length, St- As can 
be seen, in this case dN/d log u decreases less rapidly. It is 
because the total length of a given filament is usually longer 
than its spatial extent. The more convoluted a filament is, 
the larger the difference between St and S is. 

To make a more quantitative test of the universality 
of dN/dlog ly, we conduct a x^-statistics. Assuming the nu- 
merical result of dN/d log v at each redshift as one indepen- 
dent realization, we calculate the one standard deviation of 
dN/d\ogv between the four realizations, aj. One scatter er- 
ror o in the measurement of universal value is calculated 
as (T = -^(tI -|- cTp where denotes the Poisson errors. 
With this error at each bin, we evaluate the reduced x^i 
= /ni (rif: the degree of freedogm) averaged over all 
size bins. The values of Xr are listed in Table |3] where the 
best-fit values of A and B can be also found for the cases 
of two different definitions of S. As can be seen, the values 
of x^ S't'e reasonably close to unity, which proves that the 
rescaled size distribution of void filaments is almost univer- 
sal. But, it is worth noting that when the filament size is 
defined as the spatial extent, the distribution is more insen- 
sitive to redshift. 
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Figure 7. Same as FigE] but the sizes of void filaments are mea- 
sured as the sum of the lengths of all edges that constitute the 
filaments, 5^- 



Table 3. The results of the reduced test for the size distribu- 
tion of void filaments for the two different cases of the filament- 
size definition. 



size 


A 


B 




spatial extent 


0.00127 


0.569 


1.71 


total length 


0.00123 


0.859 


2.05 



5 DISCUSSION 

It is generally thought that the large scale structure of the 
universe provides a window on the initial condition of the 
early universe. In previous cosmological studies, it was the 
halo mass function that has been most highlighted as a 
statistical tool for probing co smolo gy with the large scale 
structure (Press & Schochtcr 1974: S heth fc TormenI Il999l : 
ISheth et al . 2001; Jenkins et al. 2001: Reed et al.lbOOa '). The 
most advantageous aspect of the halo mass function is that 
it can be written in a simple uni versal form, independent of 
redshift (jSheth fc TormenI 1 19991 ). Yet, the halo mass func- 
tion has a generic weakness as a cosmological probe: its de- 
pendence on the cosmological parameters comes indirectly 
from the dependence of the halo mass on the rms fiuctua- 
tion, (j{R), of the linear density field smoothed on a scale 
radius of R. 

It is desirable to develop another statistical tool 
which can overcome the weakness of the halo mass 
function. Recently, several authors have noted the 
possibilit y of probing cosm ology with cosmic voids 
(e.g.. S heth fc van de Weygaert ^04; jPark fc Leel l2007al : 
[platen et al.ll2008l ). Yet, no new statistical tool that can 
really compete with the halo mass function h as been sug- 
gested so far. For instance, in our previous work l|Park fc Led 
l2007al ). we have proposed that the cosmic voids provide an- 
other window on cosmology. Given that the void ellipticities 
are induced by the primordial tidal effect which depends on 
the background cosmology, we have shown analytically that 
the void ellipticity distribution can be used to constrain the 
cosmological parameters. 
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However, it has turned out that the void elhpticity dis- 
tribution suffers from the same weakness that the halo mass 
function has. Its dependence on cosmology is secondary, re- 
sulting from the dependence of the void scale on the linear 
rms density fluctuation, a{R). Furthermore, the void elhp- 
ticity distribution cannot be written in a universal form un- 
like the halo mass function. Therefore, true as it is that the 
void elhpticity distribution provides another way to deter- 
mine the cosmological parameters with the large scale struc- 
ture, it is not competitively compared with the halo mass 
function. 

In this work, we have for the first time quantified the 
effect of the tidal field on cosmic voids by the sizes of the 
void filaments. It is shown that the size distribution of void 
filaments is almost independent of redshift, having a simple 
universal form. A crucial implication of our result is that 
the number density of long filaments in voids should depend 
very sensitively on the background cosmology since the sizes 
of the void filaments reflect the primordial tidal field. 

Our result that the size distribution of void filaments is 
almost independent of redshift in spite of the well known fact 
that the sizes of voids grow very strongly with redshift can 
be understood as follows. It is true that the comoving sizes 
of voids keep growing as the voids expand faster than the 
rest of universe due to their extreme low density. However, 
the comoving sizes of void filaments should not necessarily 
keep growing with redshift. 

But for the effect of the tidal field, the spatial distribu- 
tion of void halos would be isotropic due to the gravitational 
rarefaction effect caused by the fast expansion of voids. In 
other words, while the gravitational rarefaction effect tends 
to make the spatial distribution of void halos isotropic, the 
tidal effect on voids tends to increase the anisotropy in the 
distribution of void halos, resulting in growth of void fila- 
ments. Therefore, as far as the tidal effect exists and coun- 
teracts the expansion effect, the sizes of void filaments will 
also grow as the sizes of their host voids grow with redshift. 

However, when the void size grows large enough to fi- 
nally reach its threshold at which the expansion effect over- 
comes the tidal effect, the degree of anisotropy in the spatial 
distribution of void halos will diminish and accordingly the 
sizes of void filaments will stop growing. Since this critical 
threshold of void size is determined by the single condition 
that the expansion effect overcomes the tidal effect, its value 
should be redshift-independent. At high-redshift when the 
sizes of voids were small, the strong dominant tidal effect 
increases the anisotropy in the spatial distribution of void 
halos, increasing the sizes of void filaments. At low-redshift 
when the sizes of voids are large, the strong expansion ef- 
fect prevents the sizes of void filaments from keep growing. 
Hence, the universal size distribution of void filaments re- 
flects the fact that the sizes of void filaments are regulated 
by the counter-balance between the expansion effect and the 
tidal effect. 

The fitting formula, eqs. ©-([S]), indicate that the num- 
ber density of void filaments, dN/dlogv, behaves like a 
power-law with power index of n = 1 in the short filament 
section while it decreases exponentially in the long filament 
section. A crucial implication of our result is that the ex- 
ponential decrease of the abundance of long-size filaments 
in voids should be very sensitive to the key cosmological 
parameters. Especially it is expected to depend sensitively 



on the amplitude of the linear power spectrum, as- In our 
previous work, we have already shown that the void elliptic- 
ity caused by the tidal effect increases with the value of erg 
l|Park fc Lee 2007a'). Accordingly, we expect that the num- 
ber density of long filaments in voids would increase as the 
value of erg increases. Unlike the void elhpticity distribu- 
tion whose dependence on erg is weak and indirect through 
its dependence on the rms density fluctuation, equation ((2)1 
suggests that the size distribution of void filaments should 
depend strongly on the power spectrum amplitude. 

To use the size distribution of void filaments as a cos- 
mological probe, however, there should be some additional 
tasks to be done in the future. First, the mass-to-light bias 
and the redshift distortion effects have to be account for. 
What one can observe is not halos in real space but galax- 
ies in redshift space. The size distribution of void filaments 
measured from the galaxy catalogs in redshift space could 
differ from the current result. Thus, it will be very necessary 
to investigate how the bias and the redshift distortion effect 
change the size distribution of void filaments. It is expected 
that the sizes of void filaments may increase in redshift space 
since the redshift space distortion effect tend to increase the 
anisotropy in the spatial distribution of void galaxies. Mean- 
while, the matter-to-light bias might decrease the sizes of 
void filaments, compensating for the redshift distortion ef- 
fect, since the massive halos are found to be distributed less 
anisotropically (Park & Lee 2009 in preparation). Hence- 
forth, the overall size distribution of void galaxy filaments 
in redshift space might be similar to that in real space, due 
to the competition between the two effects (in private com- 
munication with van de Weygaert). 

Second, a general analytic formula for the size distribu- 
tion of void filaments in an arbitrary cosmology is desirable 
to derive from physical principles. If an analytic model is 
found, it would allow us to understand the true physical 
meaning of eqs. (Hjl-Q. Third, it should be worth examin- 
ing whether the results depend on the choice of the filament- 
finding and the void-finding algorithms. As mentioned in §1, 
it has been found that different void-finders provide different 
void properties and the HV02 algorithm that we used here 
has a tenden cy to find larger vo ids than most of the other 
void-finders l|Colberg et al.l 120081 ). Furthermore, our results 
are not general ones but valid only for the specific cosmology 
of the Millennium Run simulation. Therefore, it is required 
to test the universality of the size distribution of void fila- 
ments with other void-finders and with other cosmologies as 
well. We plan to address these issues and report the results 
elsewhere in the near future. 
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